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We construct a two-band lattice model whose bands can carry the Chern numbers C — 0, ±1, ±2. 
By means of numerical exact diagonalization, we show that the most favorable situation that selects 
fractional Chern insulators (FCIs) is not necessarily the one that mimics Landau levels, namely a 
flat band with Chern number one. First, we find that the gap, measured in units of the on-site 
electron-electron repulsion, can increase by almost two orders of magnitude when the bands are 
flat and carry a Chern number C = 2 instead of C = 1. Second, we show that giving a width to 
the bands can help to stabilize a FCI. Finally, we put forward a tool to characterize the real-space 
density profile of the ground state that is useful to distinguish FCI from other competing phases of 
matter supporting charge density waves or phase separation. 



I. INTRODUCTION 

A large effort in condensed matter physics is dedi- 
cated to the classification of distinct phases of matter, 
exploring their properties, and establishing their robust- 
ness. Many phases of matter can be classified by the 
pattern in which a ground state manifold breaks sponta- 
neously and locally a symmetry of the underlying many- 
body Hamiltonian. The experimental discovery in Ref. i 
of the fractional quantum Hall effect (FQHE) in high- 
quality GaAs-A^Gai-a^As heterostructures subjected to 
very large magnetic field in 1982 was revolutionary in 
that the FQHE is associated to a family of phases of mat- 
ter that are not characterized by some pattern of sponta- 
neous symmetry breaking with a local order parameter, 
but instead by the notion of topological order. ^ ^ By anal- 
ogy with spontaneous symmetry breaking and its connec- 
tion to (Lie) group theory, a classification of topological 
order for many-body Hamiltonian has been undertaken.^ 

The study of all possible patterns of spontaneous sym- 
metry breaking or topological order, although immensely 
useful, fails to address a very important question. What 
is the mechanism by which a given pattern associated to 
a phase is established? This is a quantitative question 
that requires deciding in the simplest case of two com- 
peting terms in a many-body Hamiltonian (free energy) 
which one is the dominant one. Answering this question, 
although essential from a practical point of view, can be 
very difficult. Detailed prescriptions of the Hamiltonian 
matter so that nonuniverality rules when calculating the 
location in the phase diagram of the boundaries separat- 
ing phases of matter. In practice, a realistic answer to 
this question requires modern computing power if at all 
possible. 

A problem related to the stability of phases is as fol- 
lows. Suppose that we are given two physically different 
realizations of the same universal phase: are the most fa- 
vorable conditions for realizing the phase in one physical 
setting the same as in the other? More specifically, this 
paper analyzes the optimal situation for stabilizing a frac- 



tional Chern insulator (FCI). It also compares and con- 
trasts it to the optimal situation for stabilizing a FQHE 
in a Landau level. 

Fractional Chern insulators (FCIs) are strongly cor- 
related phases of matter for interacting fermions and 
bosons that have been found in Bloch bands with van- 
ishingor small bandwidth and nonvanishing Chern num- 
j^g^pHMl Perjxiionic FCIs share the universal properties 
of the FQHE that occurs upon partial filling of Landau 
levels. It was thus not suprising that the efforts to es- 
tablish that FCIs can be realized in model Hamiltonians 
began by mimicking as closely as possible the conditions 
for obtaining the FQHE in Landau levels. Bloch bands 
were chosen so as to try to mirror the properties of the 
Landau levels, which are special for (I) they support a 
topological attribute called a Chern number that equals 
one, and (II) they are independent of their momentum 
quantum number, i.e., they can be thought of as flat 
bands. Bloch bands can be made to satisfy properties 
(I) and (II) .13 However, even if the Bloch bands are made 
to satisfy conditions (I) and (II), it is not a priori ob- 
vious that energetics allow the formation of a fractional 
quantum Hall ground state. For example, the filling frac- 
tion 1/3 is commensurate to the lattice, hence one might 
expect that, instead of a fractional Chern state, interac- 
tions select a charge density wave as the ground state. 
Nevertheless, it has become apparent from exact diago- 
nalization studies of FCIs that the universal properties of 
the FQHE can be stabilized by starting from filling par- 
tially a (sufficiently) flat band of a Chern band insulator 
for some short-range interaction. 

That one can deviate slightly from condition (II) with- 
out destroying the FCI phase is expected, for the topo- 
logical phase is gapped and therefore stable to small per- 
turbations. In contrast, one cannot deviate continuously 
from condition (I), for the Chern number is an integer. 
The work presented below aims to analyze cases when de- 
parting from the conditions (I) and (II), associated to the 
FQHE in Landau bands, can lead to enhanced stability 
of the FCIs. 

We will present a lattice model for which the Chern 
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number can take values C = 0,±1,±2. We will show 
that it is more difficult to stabilize a FCI with C = 1 
than a FCI with C = 2 due to energetics for this model. 
We find that the gaps for a FCI at the filling fraction 
V = 1/b when C = 2 are almost two orders of magnitude 
larger than those for the most stable Chern fractional 
state that we found for the filling fraction u = 1/3 when 
C = l. 

To establish that the state at = 1/5 is a FCI and 
not some competing state that breaks spontaneously and 
locally some space-group symmetry of the lattice, we 
present in this paper a useful tool to distinguish frac- 
tional Chern states from charge density wave (CDW) 
states even when the lattices studied are rather small. 
We employ this diagnostic to map out the phase dia- 
gram that should emerge in the thermodynamic limit as 
a function of the strength of interactions, and to identify 
the nature of the phases that compete with the u = 1/5 
FCI at C = 2. In doing so, we substantiate our claim that 
a FCI can be stabilized by moving away from property 
(I) obeyed by the Landau levels. 

We then proceed to study the effects of relaxing the 
flatness condition (II). In order to probe the phase dia- 
gram of the system as a function of two parameters, one 
of which controls the size of the bandwidth, we use an- 
other numerical tool to help trace out the phase bound- 
aries of topological states using exact diagonalizations 
studies of small lattices. We look at the ratio between 
the observed gap A and the spread in energy 5 of those 
states that belong to the manifold of states that are to 
become degenerate in the thermodynamic limit. If L 
is the characteristic linear size of the lattice, this ratio 
A/S ~ exp(L/^) entails information on the correlation 
length ^ of the system. The correlation length ^ is a 
useful mean to measure the distance to the phase bound- 
ary. By combining the parametric dependence of A along 
contours at constant A/S^ we can estimate the region of 
stability of the FCI at u = 1/5 and C = 2 from exact 
diagonalization. We find that moving away from condi- 
tion (II) that is obeyed by Landau levels, i.e., including 
a non-zero bandwidth, can (up to a limit) help increase 
(and not decrease) the range of phase stability. 



II. THE MODEL 

Consider a noninteracting tight-binding model for 
fermions on the two-dimensional square lattice A made of 
Li X L2 sites and spanned by the orthonormal primitive 
lattice vectors and We impose periodic bound- 
ary conditions. The fermions have two internal degrees 
of freedom per site, which we denote as a spin degree of 
freedom s =t, |. Each hopping process is associated with 
a spin flip. The resulting noninteracting Bloch Hamilto- 
nian supports two bands and reads in momentum space 



H,:= ciB^-cTC^, (2.1a) 
Bh-i + i^fc;2 •= ^(sin + isin k^), (2.1b) 

^fc;3 *= ^1 COS ki + COS k2 + 

+ [cos(/c2 + k2) + cos(ki — k2)\ , 

(2.1c) 

where cj[, = (cj[, ^, cj[, ^) and cj[, ^ creates a fermion at mo- 
mentum k in the Brillouin zone (BZ) with spin s =t, i 
while a = (cr^, (72, 0-3) are the three Pauli matrices acting 
on spin space. The parameters t and /i^, /i = 1, • • • ,4 
are real. 

The role of the hopping parameters t and h^^ ji = 
1, • • • ,4 can be illustrated by fixing t > and expand- 
ing the Hamiltonian around the four inversion-symmetric 
momenta in the BZ k^^^^ = 7r(z,j), z,j = 0,1. To lin- 
ear order in the deviation p*^*-^^ = k — k^'^^\ i^j = 0, 1, 
from each of these four momenta, the Hamiltonian takes 
a Dirac form with masses given by 

^(ij) ^ {-lyh^ + {-iyh2 + /^s + {-ly^^h^. (2.2a) 

The model thus allows to independently control the sign 
and magnitude of the mass at each of the four Dirac 
points via the parameters /i^, /i = 1, • • • , 4. If all the 
Dirac points have a non- vanishing gap, the Chern number 
of each of the two bands is well defined. Each Dirac 
point contributes ±1/2 to the Chern number, with the 

sign depending on the chirality exp(i/c^*"^^ + 1^2*"^^) of the 
kinetic part of the Dirac Hamiltonian and the sign of the 
mass gap. The total Chern number of the lower band is 
then given by 

C=^ (-l)'^^'sgnm(^^'\ (2.2b) 

and therefore can assume the values C = ±2, ±1, in our 
two-band model (see Fig. [T]). 

In what follows, we study this model in presence 
of density-density interactions between fermions on the 
same lattice site {U) and between fermions on neighbor- 
ing sites (V) 

^int •= ^ ^ ^ Pr,s Pr,s' + ^ ^ ^ Pr,s Pr' ,s' 

(2.3) 

upon partial filling the lower band. Here, (•,•) denotes 
nearest-neighbor lattice sites and ^ := c]l ^ 5 is the 
density of fermions with spin s =t, | at site r G A. Both 
U and V are measured in units of t throughout the paper. 
We project this interaction to the Fock space built out 
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FIG. 1: (Color online) (a) Any one of the two Chern 
numbers of the two bands of the noninteracting Hamilto- 
nian ( |2.1| ) depends on the direction of the four-component 
mass vector h = (h^). When restricted to the form 
h = \h\ (cos sin ip cos 0, sin ip sin ^/v^, sin (p sin 0/V2), the 
Chern number phase diagram is obtained as shown (white: 
\C\ = 0, light: \C\ = 1, dark: \C\ = 2). b) Band struc- 



ture of the noninteracting Hamiltonian (2.1 ) with mass terms 
ij = 0, 1, defined in Eq. ( |2.2at . 
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(2.4b) 



where Uj^ ^ is the eigenvector of the lower band of the 
2x2 Bloch Hamiltonian Bj^ a in (2.1) at momentum 

while Xk second quantized operator that creates 

the corresponding state in the lower band. 

We are interested in the situation where the following 
hierarchy of energy scales applies 



(2.5) 



with W and m being the bandwidth and the band gap of 
the noninteracting band structure, respectively. The fact 
that m is the largest energy scale is taken into account 
by constraining the single-particle Hilbert space to the 
subspace spanned by the states in the lower band. Fur- 
thermore, we can impose that W is the smallest energy 
scale by a spectral flatteni ng o f the lower band of the 
single-particle Hamiltonian ( |2.1[ )^ 



^flat 



E 

keBZ 



~^k-> 



(2.6) 



which in effect projects all states in the lower band to the 
energy —1 and all states in the upper band to energy +1. 
This deformation of the Hamiltonian induces arbitrary- 
range hopping amplitudes in position space, but pre- 
serves its locality in the sense that these hopping ampli- 
tudes decay exponentially with distance. Interpolating 



between Hq and i^Q^*, via 



ffo(A) := (1 - A)F«-' + \H, 



(2.7) 



Yi 



FIG. 2: (Color online) (a) Energy eigenvalues of Hamiltonian 
( |2.8[ ) at A = measured relative to the ground state energy 
at the filling fraction z/ = l/3forL = 6x4 and iV = 8 par- 
ticles. The noninteracting parameters of Hamiltonian (2.8) 
are /i4 = and hi — h2 — —h^ — — t, corresponding to a 
noninteracting band with C — 1. The interacting parame- 
ters of Hamiltonian (2.8) are U — 1 and V — 0. K 3-fold 



quasi-degenerate ground state is observed in agreement with 
Ref 10 (b) Spectral flow induced by a flux insertion in the 
7i direction for N = 8 particles. The 3 lowest lying states 
do not mix with the would-be continuum of excited states 
and return to the same configuration after three flux quanta 
have been inserted (only the insertion of one flux quantum is 
shown). For comparison with Fig. (|3]), the FCI at u = 1/3 
has disappeared when V = 1. 



with < A < 1, will allow us to choose the bandwidth at 
will. 

We will begin with the study of the case A = 0, using 
the Hamiltonian 



H{X) :=H,{X) 



^int 



(2.8) 



by carrying out exact diagonalization studies of small sys- 
tems in Secs. |III| and |IV| In Section[V|we will study how a 
finite bandwidth (obtained by varying A) affects the frac- 
tional Chern states. We always restrict the single-particle 
Hilbert space to that of the lower band in our calcula- 
tions, i. e., w e assume the limit W/m <C U/m^V/m 
in Eq. Kb). 
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III. C = 1 AND C = 2 FRACTIONAL CHERN 
STATES 

The topological properties of the anticipated FCIs are 
believed to be described by a hierarchical Chern-Simons 
theor\E^H^ with the action 



(3.1) 

in the continuum limit. Here, a]^ a|, and are the 
two spatial components and the temporal component of 
i = 1, • • • , flavours of statistical gauge fields, A^, /i = 
0,1,2, is the electromagnetic gauge field, K is a. sym- 
metric X matrix with integer entries, Q is the 
component charge vector with integer entries, and the 
summation over the repeated indices is implied. The 
topological ground state degeneracy on the torus and the 
Hall conductivity are given by deti^ and 



(3^2) 



respectively. On the other hand, we may naively expect 
the Hall conductivity of the fractional Chern states to be 



where u is the filling fraction of the lower band. When 
the band is fully filled {u = 1), indeed the expectation 
is fulfilled. That this formula holds for generic filling 
is always permitted but it is not guaranteed. When the 
Berry curvature is not constant in the Brillouin zone, ctj^ 
does not need to be tied to u anymore.!^ 

It is known that topologically degenerate ground states 
on the torus can be transferred into one another by adi- 
abatically inserting flux through the torus After an 
insertion of an integer number of flux quanta equal to 
the degeneracy, one must recover the same state. For 
a many-body state |^) with N particles inserting a flux 
271^- in the direction is equivalent to imposing the 
twisted boundary conditions 



V ■ 



5 ' N 



J = 1, ■ ■ 



,N. 
(3.4) 



With this background to fraction al C hern states in 
mind, we proceed to discuss model ( |2.1| ). The princi- 
pal advantage of this model is that it allows to change 
the Chern number of the band from C = —2 to 2, in 
steps of 1. We first focus on the case C = 1, which can 
be realized by choosing = h2 = —h^ = —t and /14 = 
and was previously studied in Ref. |lOl In this case, the 
model hosts a 1/3 fractional Chern insulator state which 
we here include for completeness. 

Figure [2|a) shows the low-energy portion of the many- 
body spectrum and the expected 3-fold quasi-degenerate 



fractional Chern ground state when U = 1 and ^ = 0. 
The spectrum is plotted as a function of center of mass 
momentum Q. This is a good quantum number, for the 
many-body Hamiltonian shares the space-group symme- 
try of the square lattice. The fact that the states fall in 
the sectors with center of mass momenta {Qi^Q2)qxa = 
{(0,0), (2,0), (4,0)} agrees with the counting rule intro- 
duced in Ref. 7, The 3-fold qu asi-d egeneracy is consistent 
with a Chern-Simons theory (3.1) for a single species of 
gauge fields and K = 3^ Q = 1. 

Figure [2|b) shows the evolution of the spectrum when 
flux is adiabatically inserted in the direction for a 6 x 4 
site lattice with U = 1, y = 0. One observes that the 
quasi-degenerate ground state manifold indeed evolves 
independently from the excited states. This observation 
is again consistent with a = 1/3 fractional Chern state. 
Moreover, we checked that this state survives when ^ 
as long as the Chern number remains C = 1. 

With the help of Eq. (2.2b) one verifies that for suffi- 
ciently large /14 > 0, the Chern number becomes C = 2. 
In fact, the simplest set of parameters that hosts a Bloch 



band with C = 2 is hi = h2 



ho 



this set of the parameters, model (2.1) can be reinter- 



and /14 ^ 0. 



For 



preted as a layered model where two square lattices are 
superimposed on each other without any hopping term 
that connects them. We are going to argue that when 
(3.3) hi = h2 = h^ = 0, /14 = 0.7t, and U = V = 1, the many- 



body Hamiltonian H{X = 0) defined in Eq. ( |2.8[ ) realizes 
a FCI at = 1/5 filling. 

First, we show in Fig. [sj^a) that for system sizes 
6x5, 5x5, and 5x4 with N = L1L2/5 particles, 
the many body eigenvalues possess a 5-fold degenerate 
ground state. The total momentum sectors Q at which 
these states should appear for a fractional Chern state 
can be calculated with a generalization of the counting 
rule of Ref. [71 For a Li x L2 lattice with N particles 
where we choose to be a multiple of five, the momen- 
tum sectors are given by 



Qi 
Q2 



f[N{Li-b)/2^mN] modLi 
I N{L2-l)/2modL2 



(3.5a) 



with m = 0, 1, 2, 3, 4. If L2 is a multiple of five the same 
rule applies upon interchanging Qi for Q2. For the square 
lattices made of 6 x 5, 5x5, and 5x4 sites, the momentum 
sectors should fall at 



Qi 
Q2 

Qi 
Q2 

Qi 
Q2 



6x5 



5x5 



5x4 



(3.5b) 



respectively. This is exactly what we find [see Fig. [sj^a)]. 

Second, Fig. |3|b) shows the flux insertion in the 63 
direction for lattice size L = 6 x 5. The five states cor- 
responding to the ground state manifold evolve indepen- 
dently from the continuum and they only recover their 
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FIG. 3^ (Color online) (a) Energy eigenvalues of Hamilto- 
nian (2.8) at A = measured relative to the ground state 
energy at the filling fraction = 1/5 for three different sys- 
tem sizes corresponding to N — 6,5, and 4 particles (green 
diamonds, red squares, and black circles, respectively). The 
noninteracting parameters of Hamiltonian (2.8) are /i4 = 0.7 
and /i^ = for i = 1, 2, 3 corresponding to C = 2. The 
interacting parameters of Hamiltonian (2.8) are U — V — 1. 



All three systems show a 5-fold quasi- degenerate ground state 
(see text for details), (b) Spectral flow induced by a flux in- 
sertion in the 72 direction for N — Q particles. The 5 lowest 
lying states do not mix with the would-be continuum of ex- 
cited states and return to the same configuration after five 
flux quanta have been inserted (only the insertion of one flux 
quantum is shown). 



original position after five flux quanta have been inserted 
in the system. 

All these observations are consistent with a fractional 



Chern state described by the Chern-Simons theory (3.1) 
for two species of gauge fields and 



K = 



3 2 
2 3 



Q = 



(3.6) 



ilT = 3, Q = 1. Note that a single species of gauge 
fields with K = b would not explain the Hall conductivity 
(Jh = 2eV(5/i). 

Although these two pieces of evidence already point 
towards a 1/5 fractional Chern state, this state could still 
be a charge density wave (CDW). To show that this is not 
the case, we will put forward in Sec. [IV| a different tool 
to distinguish liquid states such as the fractional Chern 
states from other competing orders, namely phases with 
broken translational invariance such as CDW or phase- 
separated phases. 



IV. A TOOL TO DISCRIMINATE BETWEEN 
TOPOLOGICAL ORDER AND SPONTANEOUS 
SYMMETRY BREAKING IN SMALL LATTICES 

Both topological order and long-range order with a lo- 
cal order parameter are emergent phenomena that occur 
in the thermodynamic limit for a given dimensionality of 
space. Most model Hamiltonians are not exactly solv- 
able in the thermodynamic limit. Approximations such 
as variational methods can be used to access the thermo- 
dynamic limit, but they are uncontrolled. Exact diago- 
nalization techniques are not biased, but they are limited 
to finite sizes. The extrapolation of exact finite-size spec- 
tra to the thermodynamic limit must be undertaken with 
great care. In particular, this extrapolation can involve 
subtle dimensional crossovers (e.g., the thin torus limit 
that will be discussed below). 

In view of the effects of finite-size corrections to the 
thermodynamic limit, it is essential to use complemen- 
tary probes for topological order in exact diagonalization 
studies to argue convincingly that a candidate FCI has 
the featureless character that is demanded from a topo- 
logical fluid and to map out its boundary in the phase dia- 
gram. The particle entanglement spectrum has been used 
in Ref. ^to identify FCIs. We are going to give another 
criterion that distinguishes a fractional Chern phase from 
a competing phase that breaks spontaneously the space- 
group symmetry of the lattice. Examples of such com- 
peting phases are charge density waves and macroscopic 
phase separation. Equipped with this diagnostic, we are 
going to map out the phase diagram as a function of 
the dimensionless parameter —2<V/U < 1 of Hamilto- 
nian ( |2.8[ ) with a flat band, i.e., A = 0. We shall denote 
this Hamiltonian by H to shorten the notation. 

Let Oj. be any local operator defined for any site r from 
the lattice A. We project this operator onto the ground 
state manifold with the help of the Hermitean operator 



lim ( 



(4.1) 



We assume that the ground state manifold is n- 
dimensional and spanned by the linearly-independent 
many-body eigenstates • • • , |^^). By assumption, 

H shares the space-group symmetry of the lattice A. 
Hence, we can always choose |^ J to be a simultaneous 
eigenstate of the momentum operator with the center of 
mass where i = 1, • • • , n. We define the n x n Her- 
mitean matrix O^. with elements 



(4.2) 



Let v'q.^ be the i-th eigenstate of the matrix O^.. Its 

eigenvalue Xq^ is independent of the lattice site r as 
a consequence of translation symmetry of the Hamilto- 
nian.^^ If the space-group symmetry of the lattice A and 
of Hamiltonian H is not to be broken spontaneously by 
the ground state manifold, then it is necessary for all 
eigenvalues Aq ^ , • • • , Aq^^ of to be equal. The spread 
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of the eigenvalues measures the degree by which the sym- 
metry associated with the operator is broken. If Aq^" 
and Xq^^ are, respectively, the minimum and maximum 
eigenvalues, the quantity to monitor is the difference 



o •= 



\ max 



\ mm 



(4.3) 



We shall demand the condition SXq holds in a suit- 
able thermodynamic limit for all local operators if topo- 
logical order is to hold. 

We now apply this analysis for the case when = , 
the local density operator being defined by 



Pr 



1 



EE 

q,k s 



iqr t 



and construct the matrix with elements 



(4.4) 



(4.5) 



A set of n maps of the local fermion density in the 
ground state manifold is obtained as follows for some 



finite lattice A. We denote with 



i = 1, 



set of orthonormal eigenvectors of at some arbitrarily 



chosen site Tq and evaluate the n real functions 



,,(0t 



,(0 



,n. 



(4.6) 



The functions ^ are density maps of the n linearly 
independent combinations of the states |^^),--- , |^^) 
selected by the set of eigenvectors v^)r^^ i = 1, • • • , n. 
These functions show the variation of the local fermion 
density in position space, i.e., the finite lattice A. The ex- 
trapolation of these density maps to the thermodynamic 
limit can be used to distinguish between a ground state 
manifold that supports a charge-density wave, a phase 
separation in position space, or is featureless as would be 
expected from a FCI. We will give two examples. Exam- 
ple 1 is the case of the cross-over of a fractional Chern 
state from a featureless liquid to a charge density wave 
as the aspect ratio of the lattice is varied. Example 2 is 
the case of a phase-separated ground state obtained with 
attractive interactions between the single-particle states 
from the fiat band. 

Fractional quantum Hall states on the torus turn 
smoothly into a charge-density wave state, if the aspect 
ratio min(L;L5 ^2)/^ ijSoHlI ^.j^jg go-called thin-torus 
limit, that also exists for fractional Chern insulators, 
the occupation numbers of the charge-density wave or- 
bitals are in direct correspondence to those for the frac- 
tional quantum Hall state. For concreteness, let us con- 
sider the u = 1/5 fractional Chern state at Chern number 
2. For the aspect ratio L-^/ L2 = 1, as is the case for a lat- 
tice of 5 X 5 sites with 5 particles, all 5 topological ground 
states have the same center of mass momentum Q = 0. 
Therefore, Q{r) is the unit matrix for all r G A and all 



functions 



1, • • • 5 are independent of r. As a 



consequence, the fractional Chern state is indeed feature- 
less in this isotropic case. Upon choosing the slightly 
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n the FIG. 4: (Color online) Density profile rir^ defined in 



Eq. (4.6) for one representative many-body ground state |^^) 
among the n linearly-independent quasi- degenerate many- 
body ground states. Note the difference in the color scale 
between the plots in (a) and (b). (a) Case of the filling 
fraction z/ = 1/5 when the representative ground states for 
5 X 6, 3 X 10, and 2 x 15 lattices would turn into a fractional 
Chern insulator in the thermodynamic limit, by which the 
two linear dimensions of the lattice are much larger than the 
correlation length. In a finite lattice, a charge- density wave 
profile becomes more pronounced when the aspect ratio mim- 
icks the thin-torus limit, by which the termodynamic limit 
is taken with one of the two linear dimensions of the lattice 
comparable or smaller than the correlation length. The in- 
set shows the evolut ion of the eigenvalues \p of the matrix 
Qr-^ij defined in Eq. (4.5) as the aspect ratio of the lattice is 
changed, (b) Case of the filling fraction v — 1/5 when the 
representative ground state for a 5 x 6 lattice would turn into 
a phase-separated ground state in the thermodynamic limit, 
by which the two linear dimensions of the lattice are much 
larger than the correlation length, as a result of attractive 
nearest-neighbor interactions. 



anisotropic lattice = 5, L2 = 6 with N = Q particles, 
the eigenvalues of Q{r) lie between 0.994 and 1.007. This 
gives rise to a small density variation in position space 
of about one percent [see Fig. [4]^ a)]. Increasing further 
the anisotropy to = 2, L2 = 15 with N = Q parti- 
cles results in a spread of the eigenvalues of Q{r) between 
0.91 and 1.05 fsee Fig.ga)]. In this case, the functions 

nr\i = 1, • • • 7 5 form the expected charge-density wave 
pattern with pronounced density minima and maxima 
along the direction, as shown in Fig. [4]^ a). 

We want to contrast these results with the density pro- 
file of the ground state that emerges when attractive in- 
stead of repulsive interactions are added to the noninter- 
acting fiat band model. We choose U = 1 and vary V 
from positive to negative values for a = 5, L2 = 6 
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FIG. 5: (C olor online) (a) Quantum phase diagram of Hamil- 
tonian (2.8) at A = for fixed noninteracting parameters cor- 
responding to two flat bands with C = ±2 as a function of the 
nearest-neighbor interaction strength V for fixed on-site in- 
teraction strength U > 0. The noninteracting parameters are 
/i4 = 0.7t and — ^ for i — 1, 2, 3. There are three phases. 
For y > the zy = 1/5 fractional Chern insulator (FCI) is 
the most stable phase. For V <0 the system phase separates 
(PS) and below a critical ^ —0.6 it turns superconducting, 
(b) Superfluid phase stiffness as defined in Eq. (4.7) for the 
two orthogonal spanning directions of the lattice. 



with A/" = 6 particles. In doing so, we encounter two 
phase transitions towards different gapped ground states 
at about F - and ^ = -0.6 [see Fig. [sfb)]. 

First, we obtain a state with macroscopic phase sepa- 
ration of the fermions and quasi-degeneracy L2 = 6. The 
strongly varying density profile in real space [see Fig. 
[4]^b)] reveals that all particles cluster in a single stripe 
and the degeneracy emerges from shifting this stripe in 
the 62 direction across the lattice. 

Second, a state with quasi-degeneracy 30 = x L2 is 
obtained with a density profile that shows a stripe similar 
to the first phase. 

To characterize the physical difference between the two 
phases stabilized by an attractive V, we have computed 
the phase stiffness (see Refs. 33,.34* and.35 ) of the ground 
state 1^) against twisting the boundary conditions away 
from periodic ones by a complex phase 7^ , i = 1 , 2 ac- 
cording to Eq. (3.4). This stiffness equals the condensate 



fraction p\ of superfluid pairing in the state |^) and it is 
given by 



1,2, 



(4.7) 



7i = 



where Eq = is the many-body ground state en- 

erg}^. We find that in the phase with small negative V, 



there is no superfluid pairing, while the phase at large 
negative ^ is a one-dimensional superconductor nucle- 
ated in the stripe of clustered particles. The phase stiff- 
ness is plotted for both directions in Fig. [sjb). The large 
increase of pi when V approaches ~ —0.6 from be- 
low suggests a phase transition. The direction for which 
the phase stiffness is the largest and diverging upon ap- 
proaching the critical value ~ —0.6 from below is 
This is the direction of the stripe along which the elec- 
tron density is uniform. Hence, it is plausible to interpret 
the phase F < 1/ as a stripe-like superconducting phase. 



V. THE EFFECT OF BAND DISPERSION 

Having established the existence of fractional Chern 
ground states out of band insulators supporting the 
Chern numbers 1 and 2, we now study the effects on 
FCIs resulting from the spectrum of the noninteract- 
ing Hamiltonian not being flat. Bringing back the 
dispersion of the Bloch bands amounts to assigning a 
momentum-dependent energy penalty for occupying the 
single-particle states created (by xL) the BZ. Increas- 
ing the bandwidth of the noninteracting Bloch bands to 
some critical value must result in a phase transition to a 
new correlated ground state, say a state that supports 
long-range order or a gapless featureless ground state 
such as a Fermi liquid. For a Laughlin state at filling 

lifting gradually the degeneracy of the Landau levels 
should gradually weaken the many-body gap as we now 
argue. After projection into the lowest Landau level, the 
only remaining energy scale is the (screened) Coulomb 
interation that delivers a many-body correlation length 
^ ^ l^B ~ \/^c/ (eB) the magnetic length]. Any 
small one-body perturbation that breaks Galilean invari- 
ance, an impurity potential or a periodic potential, brings 
about a characteristic length £. The effect of this one- 
body perturbation on the distance to a gap-closing phase 
transition should depend solely on the dimensionless ra- 
tio ^/^, in which case a monotonically decreasing depen- 
dence is to be expected. Indeed, the Berry curvature of 
a Landau level is constant in momentum space. It does 
not favor any particular finite length scale. Thus, the dis- 
tance to a gap-closing phase transition upon increasing £ 
should depend solely on the dimensionless ratio In 
contrast, a FCI with flat Bloch bands has a Berry cur- 
vature that must necessarily vary in momentum space. 
The Berry curvature might thus favor certain character- 
istic lengths. The distance to a gap-closing phase transi- 
tion upon increasing i should depend on more than one 
parameter, in which case non-monotonic dependence on 
£ becomes a possibility. We are going to give supporting 
numerical evidences for a scenario by which switching on 
a finite bandwidth according to Eq. (2.7) can enhance 
the stability of a FCI. 

First, we need a measure for the stability of a phase 
that we identify by a set of quasi-degenerate ground 
states of small lattices. 
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FIG. 6: (Color online) Stability of the FCI ph ase a t = 1/5 for C = 2 as the bandwidth of the noninteracting Hamiltonian 
is changed via the parameter A defined in Eq. (2.7). (a) and (b) show the gap A between the 5-th and 6-th lowest energy 



eigenvalue of the many-body Hamiltonian and the quotient A/(5, where 6 is the spread in energy of the five lowest eigenvalues, 
respectively. The region colored blue in (a) and bright (b) is interpreted as a FCI with = 1/5. In panel (c) the evolution of 
the lower portion of the many-body spectrum with increasing A at constant /i4 = 1.05 is plotted. For all figures the number of 
particles is = 6 particles, while Li = 6, L2 = 5, /i^ = 0, i = 1, 2, 3, and [/ = F = 1.5. 



One measure for the stability of the phase is simply the 
size of the gap A above the ground state manifold, i.e., 
the difference in energy between the highest of a tuple 
of quasi-degenerate ground state energies and the lowest 
energy eigenvalue above it. 

A second measure for the stability of the phase is the 
correlation length ^ characterizing the exponential decay 
of the ground state expectation values of products of a 
pair of local operators separated by the distance r. The 
shorter ^ is, the farther is a ground state from a quantum 
phase transition. If a manifold of ground states consists 
of topologically degenerate states as is the case for the 
FCI, they become exactly degenerate in the thermody- 
namic limit with the property that no local operator can 
transform one of the ground states into another. For a 
finite system, in contrast, a splitting S in energy between 
the highest and the lowest of the quasi-degenerate ground 
states is to be generically expected. This splitting can be 
used as a measure of ^ through the Ansatz 



oc e 



(5.1) 



where L is the characteristic linear size of the system. 

For concreteness, we focus on the FCI at = 1/5 with 
a 5-fold quasi-degenerate ground state. We interpolate 
between the flat band and the original noninteracting 
Hamiltonian with the help of the family of noninteracting 
Hamiltonians (2.7) parametrized by A G [0, 1] and calcu- 
late how the stability of the candidate FCI phase changes, 
as measured by A and A/ 5. Varying A from to 1 allows 
to change the bandwidth W of the lower band relative 
to the energy scale of the interaction. The bandwidth of 
the Bloch band is thus given by ly = A Wq^ where Wq is 
the bandwidth of the noninteracting band, that depends 
on the parameters of the model, in our case. We do, 
however, project the Hilbert space to the one spanned by 
the single-particle states of the lower band for all values 
of A. 



We present in Fig. [g] both A and A/ 5 as a function 
of the bandwidth-parameter A a nd th e parameter of 
the noninteracting Hamiltonian ( |2.1| ), where the band- 
width is W = Amax{V^, \ y/2 — h^\}. While the largest 
absolute gap A is indeed obtained for the limit of flat 
bands A = in this parameter space, we observe that 
for a large range of values of both the gap A as well 
as the quotient A/S can be substantially increased as A 
becomes larger. We interpret this as an increase in the 
stability of the FCI phase with increasing bandwidth of 
the noninteracting band and exemplify this scenario for 
/14 = 1.05 in Fig. ^c). Furthermore, we observe that 
the FCI is stable against a substantial bandwidth of the 
order of the interaction energy scale for ~ 0.7. 



VI. CONCLUSIONS 

The search of materials that realize FCIs requires the 
hierarchy of energy scales 



ly < V < m. 



(6.1) 



with the largest and smallest energy scale the band gap m 
and the bandwidth W of the noninteracting band struc- 
ture, respectively, while the intermediary energy scale V 
is the characteristic interaction energy. This work as- 
sumes the limit 



W/m < V/m 



(6.2) 



and deals with the energetic question of what is the op- 
timum ratio of W/V needed to stabilize FCIs. 

In this work we have established that fractional Chern 
insulator phases can be stabilized in situations that differ 
considerably from the ordinary Landau level paradigm, 
namely a flat band with Chern number one. We have 
characterized the phase diagram of a model that hosts 
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Bloch bands with Chern number either one or two by 
means of exact diagonahzation. Together with the usual 
characterizations via flux insertion and the counting of 
degenerate ground states, we have mapped out the real 
space density profile of several phases that appear in this 
model at filling = 1/5 for C = 2: charge density wave, 
phase separation, and fractional Chern insulator. We 
showed that the gaps for the = 1/5 fractional Chern 
insulating state with C = 2 are more than one order of 
magnitude larger than those for the most stable state 
that we found in the model for C = 1, at u = 1/3. 
From these studies we conclude that the most favorable 
conditions to stabilize a FCI need not be those with C = 
1. 

By analyzing the phase boundaries of the z/ = 1/5 frac- 
tional Chern insulating state for C = 2, we have found 
instances for which the flat band condition is not opti- 
mal. Instead a nonvanishing bandwidth is. In these ex- 
amples, the stability of the FCI first increases with band- 
width, reaches a maximum, and then decreases again as 
too large a bandwidth compared to the interaction energy 
scale disfavors the topological state. We have thus shown 
that it is possible to enhance the stability of a fractional 
Chern insulator by moving away from the condition of 
fiat bands. 

The energetic question of when electron-electron inter- 



actions can select strongly correlated phases of matter 
supporting topological order is much more subtle when 
electrons populate the Bloch levels of a Chern band insu- 
lator instead of the lowest Landau level. This difference 
manifests itself with the fact that the Berry curvature is 
never constant throughout the Brillouin zone in the for- 
mer case, whereas it is perfectly flat for the lowest Lan- 
dau level in the latter case. In this sense, the competition 
between the kinetic energy and the electron-electron in- 
teraction is much richer for FCIs than it is for the FQHE 
from the lowest Landau level. These subtleties highlight 
the importance of the interplay between topology and 
energetics. 
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